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It has been proposed that the superconductivity of Sr2Ru04 is characterized by pairing that is 
unconventional and, furthermore, spontaneously breaks time-reversal symmetry. However, one of 
the key expected consequences, viz., that the ground state should exhibit chiral charge currents 
localized near the boundaries of the sample, has not been observed, to date. We explore an alter- 
native implication of time-reversal symmetry breaking: the existence of walls between domains of 
opposing chirality. Via a general phenomenological approach, we derive an effective description of 
the superconductivity in terms of the relevant topological variables (i.e., domain walls and vortices). 
Hence, by specializing to the in the in-plane rotationally invariant limit, we show that a domain 
wall that is translationally invariant along the z axis and includes a bend through an angle is 
accompanied by a nonintegral (and possibly nonquantized) magnetic bend flux of ((S/tt) + n)$o, 
with integral n, that penetrates the superconductor, localized near the bend. We generalize this 
result to the situation in which gauge transformations and rotations about the z axis are degenerate 
transformations of the chiral superconducting order. On the basis of the specialized result and its 
generalization, we note that any observation of localized, nonquantized flux penetrating a 2-axis sur- 
face (e.g., via scanned-probe magnetic imaging) can be interpreted in terms of the presence of bent 
walls between domains of opposing chirality, and hence is suggestive of the existence of time-reversal 
symmetry-breaking superconductivity. 



I. INTRODUCTION 

Recently, there has been developing excitement regard- 
ing the nature of the superconducting state of the crys- 
talline compound Sr2Ru04 1 . This is because, as in 
superfluid ^He, the superconductivity has been proposed 
to be unconventional, having Cooper pairs of the triplet 
type PJ-H] ■ In addition, recently obtained evidence for the 
existence of half-quantum vortex structures [Si , which are 
expected to support zero-energy Majorana modes [SI [7], 
suggests that Sr2Ru04 could potentially be used as a 
host medium for topological quantum computing [3] ■ 
However, questions remain concerning the structure of 
the pairing state (see, e.g., Refs. [TOl [TTj) and, in par- 
ticular, whether the superconductivity does indeed spon- 
taneously break time-reversal symmetry, and would thus 
form a chiral state [I2l - [T4] . In particular, the theoreti- 
cal prediction (see e.g. [T3]) that the ground state should 
exhibit chiral charge-currents localized near the bound- 
aries of the sample, has not been verified experimentally, 
to date, despite considerable efforts [16-18^. Moreover, if 
time-reversal symmetry were broken by the superconduc- 
tivity of Sr2Ru04 then — in addition to vortices — domain 
walls that separate regions of opposing chirality would 
enter as a new topological feature of the theory [T^] . 

For conventional superconductivity, the phenomeno- 
logical approaches of London ^20^ and of Ginzburg and 
Landau [H] predated the formulation of the microscopic 
theory, due to Bardeen, Cooper, and Schrieffer [22]. 
For unconventional superconductors, including those in 
which time-reversal symmetry is spontaneously broken, 
it is likewise possible to make progress phenomenologi- 



cally, without invoking detailed information about any 
specific microscopies [THl HSl [H] . This is the approach 
that we adopt in the present paper, as we explore certain 
specific features of time-reversal symmetry-broken states: 
(i) the possibility that there are walls between domains 
of opposing chirality, (ii) the threading of these walls by 
magnetic flux, and (iii) the fact that this flux may pen- 
etrate in nonquantized amounts [251 I26j . Lack of flux 
quantization has been discussed in related settings, such 
as superfluid condensates of ionized hydrogen [37] , as well 
as time-reversal symmetry-broken superconductors that 
feature spin-polarization j2Sj, disclinations [21], or inter- 
secting grain-boundaries between crystallites [30] . 

Our central result is as follows: nonintegral (and even 
nonquantized) multiples of the superconducting quan- 
tum of magnetic flux penetrate time-reversal symmetry- 
breaking superconductors, localized near bends in walls 
between chiral domains. We first obtain this result via 
an effective description in terms of domain walls and vor- 
tices, which shows that (in the special case of the crys- 
tallographically in-plane rotationally invariant limit) a 
domain wall that is translationally invariant along the z 
axis and bends through an angle is accompanied by a 
net flux (which we term "bend flux") of ((O/tt) -f ri)$Oj 
localized in the bend region, for arbitrary integral n. We 
then generalize this result to the situation in which gauge 
transformations and rotations about the z axis are degen- 
erate transformations of the maximally chiral supercon- 
ducting order (i.e., are transformations that have equiva- 
lent impacts) . If the rotational symmetry is broken down 
to discrete tetragonal symmetry, our central result re- 
mains valid for the particular case of a domain wall bent 
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through 7r/2 radians. 

Our speciaUzed and more general results indicate that 
observations of localized, nonquantized flux penetrating 
a z-axis surface (e.g., via scanned-probe magnetic imag- 
ing) , could potentially be interpreted in terms of the pres- 
ence of bent walls separating domains of opposing chiral- 
ity, and hence would be suggestive of the existence of 
time-reversal symmetry-breaking superconductivity. Al- 
ternatively, if localized nonquantized flux is not observed 
to penetrate a z-axis surface, this would suggest that ei- 
ther (i) domain walls are not present, or (ii) domain walls 
are present, but are arranged in a parallel array and thus 
are not bent. 

This paper is organized as follows. In Section |ll] we 
review the structure of the Ginzburg-Landau order pa- 
rameter appropriate for unconventional superconductiv- 
ity with broken time-reversal symmetry, along with the 
corresponding Ginzburg-Landau free-energy functional. 
In Section |III| we analyze this free energy via an ex- 
tension of the London limit, in which we exchange the 
Ginzburg-Landau order-parameter description for a re- 
duced description in terms of the collection of spatially 
varying "phase-like" fields that for homogeneous con- 
figurations would parametrize the space of equilibrium 
states. The extension amounts to taking the limit in 
which the domain walls are vanishingly thin, compared 
with the London penetration depth. (It should be noted 
that Heeb and Agterberg, in Ref. [ST], attach a differ- 
ent meaning to the term "extended London limit'.') Tak- 
ing this limit enables us to focus on the structure and 
implications of the topological excitations of the order 
parameter, which have the form of vortices and domain 
walls, and to relate the densities of these excitations to 
singularities in the phase-like fields. In Section |IV| we 
return to the free energy, expressing it in terms of these 
excitation densities, and in Section [V] we use this frame- 
work to determine the spatial distribution of magnetiza- 
tion associated with domain-wall topological excitations, 
specifically walls that contain bends. In this section we 
also show that the thin-domain-wall limit is not essen- 
tial, in the sense that the key features of our results, 
such as the threading of bent domain walls by nonquan- 
tized amounts of magnetic flux, continue to hold, even 
when this limit is relaxed. In Section IVII we consider 
three experimental settings in which nonquantized flux 
may be observed; positive results in any one of them 
would provide evidence for the existence of time-reversal 
symmetry-breaking superconductivity. Finally, in Sec- 
tion |VII| we summarize our key results and their impli- 
cations. Some technical details are relegated to a pair of 
appendices. 



II. PHENOMENOLOGICAL THEORY OF 
UNCONVENTIONAL SUPERCONDUCTIVITY 

In this section we describe the phenomenological the- 
ory of superconductivity, on which our analysis is rooted. 



This approach is based on the notion of an appropri- 
ate superconducting order parameter, along with general 
symmetry considerations, and thus can be explored in- 
dependently of any specific microscopic details. The or- 
der parameter transforms under the full symmetry group 
of the physical system, and thus provides a representa- 
tion of this symmetry group. (It is common, in the con- 
text of planar superconductors, for the point-group as- 
pect of this symmetry to be the tetragonal group D4h, 
reflecting the underlying electronic and atomic struc- 
ture of the crystalline material.) For Sr2Ru04, the ma- 
terial on which we shall focus, it is known that, for a 
range of temperatures close to the superconducting tran- 
sition temperature Tc, the superconducting properties 
are nearly isotropic with respect to rotations about the 
z axis |32l I33j (i.e., the direction perpendicular to the 
Ru02 planes), and is only weakly tetragonal about this 
axis. Although, accordingly, the initial focus of our anal- 
ysis will be on the isotropic limit (which we term the 
in-plane rotationally or S0(2)^-invariant limit), we do 
subsequently address the cases in which the symmetry is 
lowered to the discrete group D4fi (and also, parentheti- 
cally, the group Dq^)- At the outset, we therefore retain 
generality by determining the representation furnished 
by the superconducting order parameter appropriate to 
D411 symmetry, motivated by the relevance of this group 
to Sr2Ru04. 

We now determine the appropriate representation of 
the superconducting order parameter, bearing in mind 
the foregoing symmetry considerations. This choice of 
representation is made according to the following three 
simplifying assumptions: (i) The ground state of the su- 
perconducting order should transform trivially under lat- 
tice translations. Thus, at the lengthscales relevant for a 
phenomenological description such as the one used here, 
the ground state of the superconducting order is transla- 
tionally invariant, (ii) The representation of the symme- 
try group should be irreducible. This is justified in the 
case of Sr2Ru04 as, in the absence of an applied magnetic 
field, only one superconducting transition seems to be ob- 
served. (Recent experiments on Sr2Ru04 under uniaxial 
pressure do, however, indicate the possibility of a sec- 
ond transition [H].) (iii) The representation should allow 
for the possibility that the superconducting state spon- 
taneously breaks time-reversal symmetry. This would re- 
quire that the dimension of the representation be greater 
than unity. Taken together, these assumptions fix the 
order parameter to transform according to the Fs rep- 
resentation [23], which is two-dimensional [35 . Accord- 
ingly, the order parameter is the complex- valued, two- 
component field ?/a(r), where the index a runs through 
the corresponding basis functions of the representation 
(i.e., X and Y) which, in general, depends on the three- 
dimensional position vector r. To simplify our analysis, 
we consider superconducting states that are translation- 
ally invariant along the z axis, thus rendering the physical 
problem effectively two-dimensional. Provided we apply 
external magnetic fields that are oriented along the z di- 
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rection (i.e., H = Hz), this is an option, owing to the 
intrinsic translational invariance of the system along the 
z direction. These requirements, taken together, then 
dictate that the appropriate Ginzburg-Landau free en- 
ergy functional governing the representation is given 
by [HIES] 



(1) 



F'W] = J d\'{K.:,,, {D'M*{D'M -All': il 

+ \B'a,caV': Vt V'c V', + ^l(V' X A') - Hf }. 

Here, two-dimensional summations are implied over the 
repeated indices a, b, c, d, and the covariant derivative is 
defined via D' := V' — 27rzA'/$o, where $0 is the super- 
conducting flux quantum hc/2e. 

In Eq. ([1]) the primed variables are dimensionful. We 
now define relevant scales of length and energy, and use 
them to introduce convenient dimensionless variables, 
which we use throughout the remainder of the paper and 
which we write without primes. As a first step, for the 
coefficient tensors IC' and B' we define the dimensionful 
scale factors /C and B, which we then use to construct 
the dimensionless tensors IC /C'//C and B := B' /B. In 
the in-plane rotationally invariant limit, symmetry con- 
siderations dictate that K, and B can be parametrized in 
the following way: 



Babcd — lachd + ^-^'^ac^^bd^ 



K-abcd = lachd + ^J-EacEbd + ^^'-^ac^^bd^ 



(2) 



where a summation from f to 2 is implied over the re- 
peated index S, the three real parameters {a, fi, r} are, 
in principle, temperature dependent, and the constant 
tensors {I,E, , M^} are defined via 



/ 



1 

1 
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1 
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The choice of tensor decomposition in Eq. ([2| is moti- 
vated by the observation that under 80(2)^ rotations 
I and E each transform trivially, whereas the pair 
{M^ , M^} mix. If the symmetry were not 80(2)^ but 
D4h, the coefficients of the the M^M^ and M'^M^ terms 
may be unequal; however, for D^h symmetry they would 
remain the same. 

We choose a natural scale for the order parameter, viz., 
770 = {A/By^^, and use it make the definition rj :— rj' /r]^. 
We then define the two lengthscales: (i) the penetra- 
tion depth A := 'i>o/(327r^77QXi)^/^, which characterizes 
the lengthscale for variations of the magnetic field; and 
(ii) the coherence length ^ := {Kl/A)^^"^, which charac- 
terizes the lengthscale for variations in the amplitude of 
the order parameter. We then scale all lengths by A, via 
(r, V,D) := (r'/A, AV, AD')- We also make the con- 
ventional definition of the Ginzburg-Landau parameter 



AC := A/^. Next, we define the dimensionless vector po- 
tential A, applied magnetic field H, and magnetic flux 
$ via (27rA747$o,27rA2H7$o,27r$'/$o)- 
We note that, with this choice of units, a flux equal to a 
fiux quantum has the dimensionless value 2tt. As a final 
step in the construction of the dimensionless variables 
we choose as a scale for the free-energy density the value 
/o = 2/CryQ/A^, using which we arrive at the dimension- 
less free energy via F :— F' /X^/q. It will be convenient 
for us to separate contributions to the dimensionless free- 
energy density into two groups: the "London type" terms 
/l and the and potential terms /p, respectively defined 



via 



h = llCabcd{DaVbnDcVd) + ^|(V X A) - Hp, (4a) 



/P = ^fi^i-VlVa + ^BabcdVl Vb Vc Vd), 



so that 



dMh + /pot}. 



(4b) 



(4c) 



One way in which F differs from the conventional 
Ginzburg-Landau free energy is that the tensors IC and B 
allow for a larger number of material-dependent param- 
eters, the latter free energy having only a single such pa- 
rameter, viz., K. It is possible to estimate these additional 
parameters under the assumptions of weak coupling and 
a cylindrical Fermi surface [35], and this results in the 
following values: (^,t, ct) = (0,1,1). However, due to 
the presence in Sr2Ru04 of effects such as multiple elec- 
tronic energy bands [S^ , spin-orbit interactions [S^ , and 
chiral currents [15] , the parameters of a Ginzburg-Landau 
theory that incorporates such effects self-consistently are 
expected to be modified from their weak-coupling val- 
ues, perhaps significantly. Thus, we shall not limit our 
analysis to the weak-coupling values of these parameters. 

In the section that follows, we analyze the potential 
terms of the phenomenological free energy, Eq. (4b), and. 



specifically, review how its structure leads to both vor- 
tices and domain walls. In particular, we derive the vor- 
tex and domain- wall densities in terms of the "phase-like" 
variables; in the subsequent section. Section [IV[ we con- 
struct the effective free energy in terms of topological 
variables, such as the vortex and domain- wall densities. 



III. TOPOLOGICAL FIELD CONFIGURATIONS 

As is well known, for many purposes, the state of an 
ordered phase can be adequately specified via an order- 
parameter field that takes values lying in the subspace of 
degenerate homogeneous equilibrium states TZ (see, e.g., 
Ref. [33] )• If) as an example, different regions of a sam- 
ple were to adopt distinct such values, it can — depending 
on the structure of the order parameter — be possible for 
the system to become trapped into order-parameter con- 
figurations that possess topologically stable defects [39) . 



These are spatially varying configurations of the order pa- 
rameter that cannot be removed via local deformations. 
The framework of homotopy groups of TZ enables one to 
identify and classify the possible topologically stable de- 
fects. 

As is also well known (see, e.g., Refs. [Il l40ti42] ). 
there can be a rich interplay between the topological 
features of the (bosonic) order-parameter fields that de- 
scribe ordered phases and the qualitative character of 
any fermionic particles moving in the presence of such 
order-parameter fields. However, in the present work we 
shall only consider the topological features of the appro- 
priate order-parameter field, leaving for future work the 
analysis of its implications for the motion of fermions. 

To determine TZ for the present problem, we follow 
the standard approach (see, e.g., Ref. [53]) and analyze 
the structure of the potential terms of Eq. ([T]) . To sim- 
plify the analysis, it is useful to parametrize the two 
complex fields of the superconducting order parameter 
*7 — {vxiVy) in terms of four real fields {|??|, 0, 7, 13} that 
transform simply under the operations of the symmetry 
group: 



cos 



(/3/2) \ 
isin(/3/2) ) 



cos 7 — sm 7 
sin 7 cos 7 



(5a) 
(5b) 



Now, elements of the product group U(l) x 80(2)^ 
of gauge transformations and z-axis rotations can be 
parametrized via a phase angle 9' and a rotation angle 7'. 
Under such elements, the order parameter r) transforms 
as 



r7(|r;|,^,7,/3)^r7(|r;|,0 + 0',7 + 7',/3); 



(6) 



under time reversal, tj transforms as 

r,(|,7|, 0, 7, /3)^r,(|,7|, -0,7,-/3). (7) 

Thus, we see that the parametrization of r), Eq. |5a[ is 
given in terms of an amplitude |?7| and phase 9 that are 
similar to those used in conventional superconductivity, 
but also two angular variables, 7 and that respec- 
tively characterize the additional nontrivial 80(2)^ and 
time-reversal structure associated with the version of un- 
conventional order under consideration (but see 1431 ) . In 
terms of the parametrization given in Eq. (5a), the po- 
tential terms ( 4b ) become 



/p = -yl^P+'^l^l'+8^H'cos2/3, (8) 

in which we have introduced the dimensionless length 
L ~ a^^/^/K, which will turn out to determine the 
domain-wall width. As required by 80(2)^ and time- 
reversal invariance, these potential terms are indepen- 
dent of 7, as well as being even functions of /3. If 
the symmetry were reduced to 04^, there would be 
the possibility of an additional term, proportional to 




FIG. 1. Visualizations of the order parameter space. In the 
standard London limit, in which k — >■ 00, the order parameter 
is restricted to the coset space of configurations x S'^ /Z2 
[see the discussion following Eq. (|8|]. When, in addition, the 
Ginzburg-Landau parameter a of Eq. ([2| is positive, the north 
and south poles of the sphere become free-energy minima, 
and are thus energetically preferred, relative to the equator. 
In the thin domain- wall limit (i.e., L — >■ 0), which we term 
the extended London limit, this preference is extreme. In this 
latter case, the coset space describing degenerate minimum- 
energy configurations becomes x — }. (color online) 



lyyl^ cos(47) cos^(/3). In the present setting, to achieve 
the standard London limit, in which the magnitude of 
the order parameter jryj is fixed at unity, we take the joint 
limit (k, (t) — (cx),0), keeping L fixed. In this limit, the 
structure of the order-parameter space can be visualized 
as being the product of (i) a circle, corresponding to the 
gauge degree of freedom 9, and (ii) a sphere, correspond- 
ing to the angular variables {7,/?}) (see Fig. [T]). The 
parameters {9, 7) and (0 -I- tt, 7 -t- tt) give identical val- 
ues of the order parameter [see Eq. ( 5a )] , and therefore 
correspond to physically identical configurations. 

As we aim to discuss states having time-reversal sym- 
metry breaking, we have assumed the Ginzburg-Landau 
parameter a (defined in Eq. [2]) to be positive as, for 
sufficiently weak 80(2)^ symmetry-breaking, this favors 
states in which f3 = ±7r/2. These states are related by 
time-reversal symmetry, and can be visualized as lying 
at the poles of the order-parameter sphere (see Fig. [I]) . 
In these states, the order parameter takes the form 



V\l3=±7:/2 = me 



V2 



1 

±i 



(9) 



If, across the entire system, the state had chirality f3 = 
7r/2, only a single, position dependent, phase field would 
be required to describe low-energy excitations away from 
equilibrium (and, similarly, if the state had only chirality 
/? = —tt/2). On the other hand, to describe low-energy 
excitations featuring both chiralities, as well as the "do- 
main walls" between them (through which /? changes be- 
tween ±7r/2), a pair of position-dependent phase fields, 
9{r) and 7(r), is required. It will often be convenient to 
exchange these fields for the pair 



= (r) :=0(r)T7(r). 



(10) 
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From Eq. ([s]) we can see that within a domain of maxi- 
mally positive (or maximally negative) /3, the free energy 
does not depend on (or 0^), and this remains true 
even after weak S0(2)z symmetry-breaking terms are 
included in /p. Consequently, the subspace of energy- 
degenerate homogeneous equilibrium states is discon- 
nected, being composed, topologically, of two circles, 
which exchange under time-reversal, i.e. TZ = S^x — } 
(see Fig. [I]) . This order-parameter space combines two 
of the most familiar order-parameter spaces: the 
of conventional superconductivity /superfluidity, and the 
{+, —} of Ising magnetism. 

In general, to analyze the topological structure of order 
parameters, we consider their homotopy groups 7r„(72.) 
associated with TZ. For the specific case of 72. = S*^ x 
{-h, — }, since each connected piece is isomorphic to 
the first homotopy group tti{TZ) ^ Z. This implies that 
a domain of a given chirality can exhibit vortex singular- 
ities, as, e.g., in the case of conventional superconductiv- 
ity. As the space TZ is topologically disconnected, the ze- 
roth homotopy group is also nontrivial, i.e., ttqITZ) = Z2; 
this implies the possibility of domain walls, which sep- 
arate regions of opposing chirality. (Domain walls are 
common features of systems in which the order parame- 
ter is discrete as in the Ising case.) The Z2 value of ttq 
indicates that domain walls annihilate one another. We 
remind the reader that order parameters for which 7r„(7?,) 
is nontrivial support topological defects of co-dimension 
n + 1. Thus, in the effectively two-dimensional (real) 
space that we are considering, vortices points and do- 
main walls are lines. 

The domain walls and vortices determine the quali- 
tative structure of order-parameter field configurations; 
e.g., vortices in a domain of positive or negative chirality 
correspond to topological singularities in ^'''(r) or 0~{r). 
In particular, when there are a total of vortices at 
positions {R^}^ having vorticities {q^}^ interior to 
the positive- (or negative-) chirality domain, the singu- 
larities of 9^ can be characterized by the local vortex 
density (scaler) fields and , which are defined via 



27rp^(r) ■=Eab'^a'^b 



27r 



q^S{r-R^). (11) 



Here and elsewhere in this paper, the Dirac delta func- 
tions S are are softened on an appropriate lengthscale; 
for vortices it is the vortex core diameter. 

Domain walls also have implications, but for the qual- 
itative structure of /3(r). In two spatial dimensions, do- 
main walls are lines, and a collection of N then can be 
characterized by specifying their trajectories {Rn{s)}n=i 
as functions of an arclength parameter — s„ < s < s„. By 
requiring, in addition, that the unit vector normal to the 
domain wall fiais), which is related to the domain wall 
trajectory via 

ha{s) = {cos (j){s), sin (j){s)) a = -EabdsRb{s), (12) 





FIG. 2. (a) A large region of positive chirality having an 
internal-island domain of negative chirality. As the positive- 
chirality region is multiply connected, the winding of S"*" 
around the contour Ci is an independent topological vari- 
able, (b) Six, connected chiral regions (including the exterior, 
positive-chirality region) . For every multiply connected chiral 
region, there is an independent topological variable associated 
with each non-contractible loop. As an example, the multi- 
ply connected negative-chirality region has two independent, 
non-contractible loops C2 and C3. Each of these loops is as- 
sociated with an independent winding number for 6~ . (color 
online) 



point from the negative towards the positive chiral do- 
main, the sense of the vector tangent to the domain wall, 
dsR{s), is determined. It is natural to associate the loca- 
tions of the domain walls with the zeros of the field /3(r); 
for a given set of domain walls, the equilibrium form of 
/3(r) interpolates smoothly — with a solitonic form whose 
thickness is then the domain wall width — between regions 
in which it is essentially uniform and equal either to 7r/2 
or to — 7r/2. Such structure can be characterized via a 
domain-wall density (vector) field pdw, defined via 



Pdw(r) -Vsin;3dw(r) 



J2 / dsh{s)S{r~Rn{s))- 

(13) 

Here, the delta function is softened on the lengthscale of 
the domain-wall width. We shall make use of the vortex 
and domain wall densities given in Eqs. ( 11|13 1 in Sec. IV 
in the construction of an effective free energy for the the 
topological variables. 

It should be recognized, however, that these densities 
do not, by themselves, fully specify the topological struc- 
ture of the order-parameter field. To see this, note, e.g., 
that a single island chiral domain affects the topology 
of the surrounding domain by rendering it multiply con- 
nected. Thus, to fully specify the topological structure 
of the 6'+(r) and 0~{y) fields — in addition to specifying 
the location and vorticity of the individual vortices that 
lie within the respective domains — the global winding of 
of 0^(r) or 0~{y) must be specified around a loop that 
encircles each island (see Fig. [2^) . To generalize to chi- 
ral domain structures that involve islands within islands, 
we note that to fully specify the topological structure 
of each positive (negative) connected chiral domain re- 
gion, one must specify the winding of the 6^ (^ ) field 
around each independent non-contractible loop of that 
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FIG. 3. Three types of domain wall loops (black lines). The 
gray line segments intersecting the domain wall indicate the 
local direction of the planar vector (cosr(s), sinr(s)) (i.e., 
the strongest-pairing direction). Each loop is characterized 
by the three winding numbers (n+,n_,n^ ) [see Eq. fTsfTel )]: 
(a) (-1,1,1); (b) (-2,0,1); and (c) (1,0,-1/2). These do- 
main wall loops are topologically equivalent to vortices with 
the following properties (a) singly quantized and rotation- 
ally invariant; (b) coreless and rotationally invariant; and 
(c) singly quantized and coreless. (color online) 



region (see Fig. |2|d) f?3] . 

In the remainder of this section we explain the connec- 
tion between the 7(r) field on a single domain-wall loop 
surrounding an island and the determination of the global 
freedom to wind possessed by the multiply-connected re- 
gion exterior to the island. We also examine various situ- 
ations involving individual chiral-domain islands, chosen 
to illustrate the physical roles played by the 7(r) field. 
Before doing this, we remark that on any domain-wall 
line R{s) (i.e., a locus of points on which /3 = 0), the 
order parameter takes the form 



cos 7 
sin 7 



(14) 



Thus we see that the field 7 evaluated along a domain 
wall trajectory -R(s), defines an angular variable T{s) := 
-f[R{s)) for each value of the arclength parameter s. We 
furthermore see that the function T(s) determines the 
structure of the order parameter along the domain wall- 
line, specifically via the planar vector (cosr(s),sinr(s)). 

To illustrate the physical role played by r(s), we begin 
by considering the special case of a simply-connected chi- 
ral domain, bordered by a domain wall, and thus interior 
to a larger surrounding region of opposing chirality. Two 
topological numbers (ne,n^), defined as follows, charac- 
terize the winding of the order parameter along paths 
that lie near to the domain wall: 



2Trne ■— f d9 , 

J dw 

<y dw 



(15) 



where "dw" indicated that the line integrals are evalu- 
ated along the domain wall, and the orientation of these 
integrals is taken to be counterclockwise. From Eq. (10 1, 



we see that {ng, n^) provide the same information as the 
two topological numbers (n+ , n_ ) , defined via 



2-nn± := i> d6^ 

I dw 



(16) 



specifically, n± = ng ^ n-y. We note that ng and n-^ are 
either both integral or both half-integral [so that t] is 
single- valued; see Eq. (14|], and thus that n+ and n_ are 
both integral. 

For the sake of definiteness, we consider a domain of 
negative chirality that constitutes a simply-connected is- 
land within a larger, positive chiral domain. The positive 
domain is then rendered multiply connected; see Figs. [2] 
and^ Each of the winding numbers {n^,n-,nj) corre- 
sponds to a distinct physical property: 

• The winding number of the exterior (positive) 
domain essentially determines, via $tot = ^0'^'+, 
the total flux $tot through an area that extends 
beyond the region bounded by the domain wall by 
a few penetration depths. 

• The winding number n_ of the interior (negative) 
domain is the net number of vortices in the interior 
domain; if n_ — then the domain-wall loop can 
be coreless, i.e., there is no topological requirement 
that there exist locations where |r;| vanish. 

• Whether or not the winding number n-y is unity 
determines whether or not the the superconduct- 
ing order can be rotationally invariant around a 
circular domain wall. 

Importantly, as we previously noted in this section, by 
specifying the vorticial content in the interior and ex- 
terior domains, the interior winding number is uniquely 
determined, whereas the exterior winding number is not. 
However, if in addition to the vorticial content, the value 
of n-y is specified, the winding of the exterior domain is 
also determined. In the particular case under considera- 
tion, — n_ — 2nry. 

Thus, plays a dual role, determining both the total 
fiux <I>tot, via its influence on n+ as well as whether or not 
the superconducting order can be rotationally invariant. 

We pause to make two remarks concerning issues of 
energetics. First, as a domain wall has finite energy-cost 
per unit length, to reduce its energy a domain wall loop 
may shrink in size. When viewed on a lengthscale much 
larger than the domain size, a small domain wall loop ap- 
pears topologically equivalent to a vortex [45 , and thus 
provides an alternative description of the various vortex 
structures that can occur in superconductors with bro- 
ken time-reversal symmetry (see, e.g., Refs. [3111461147] ). 
Second, in the limit in which the free energy is rotation- 
ally invariant and k is large, energy considerations prefer 
domain- wall loops that are singly quantized (i.e., contain 
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flux $o)j rotationally invariant, and coreless. However, 
owing to the linear relationship between n+, n_, and n^, 
these preferences cannot all be simultaneously satisfied. 
Compromise order-parameter configurations result from 
this frustration; we show in Fig. [3] examples of the three 
types of vortices that satisfy two of the three preferences. 
Which particular type of vortex is preferred, energeti- 
cally, will depend on the details of the parameters in the 
Ginzburg-Landau theory. 

To illustrate this frustration and the dual physical role 
played by n-^, we now consider two of the three small 
domain-wall loops that are favored energetically in the 
extreme London limit. In this limit, energetic consid- 
erations allow only coreless vortices, e.g., n_ — 0, and 
this implies that $tot — — 2$o'^7- Thus we see that, in 
the extreme London limit, if we also impose rotational 
invariance, namely = 1, we effect the magnetic prop- 
erties of the vortex, requiring the vortex to be doubly 
quantized (i.e., contain 2(f>o of flux) and flxing the sense 
of the magnetic flux. Conversely, if we fix the vortex to 
be singly quantized then, in the extreme London limit, 
the vortex would not be rotationally invariant. This in- 
terplay between the rotational and magnetic structure, 
perpetrated by the dual nature of the 7 field, underpins 
the central results of this work. 

So far, we have established that, when taken together 
with vorticial content, fixes the overall winding of the 
order parameter along a non-contractible loop within a 
multiply-connected chiral domain. However, to describe 
the local structure of a domain wall [see Eq. (14|], it 
is necessary to specify the local value of 7 [viz. r(s)] 
along the domain wall. As we shall see in Sees. |IV| and 
VB the local behavior of T{s) also plays a role in de- 



termining the flux that penetrates through the domain 
wall locally [551 US]- Thus, in order to develop a local 
description of the superconductor, the natural degree of 
freedom to use — for specifying the additional topological 
structure afforded by the presence of multiply-connected 
regions — is r(s) rather than n^. 

In the following section, by starting with the Ginzburg- 
Landau free energy, we construct an effective local free 
energy in terms of the topological variables. Specifically, 
we show that, in addition to the vortex density and do- 
main wall density, the free energy depends on a third 
topological variables, viz., the value of r(s) along domain 
walls. 



IV. EFFECTIVE FREE ENERGY IN TERMS OF 
TOPOLOGICAL DESCRIPTORS AND THE 
EXTENDED LONDON LIMIT 



tion, the focus is on the dependence of the free energy 
on the topological variables [i.e., the locations of the vor- 
tices and domain wall, as well as r(s)]. The presumption 
is that the degrees of freedom associated with exciting 
the order parameter around the state of minimum free 
energy within a fixed topological sector [defined via the 
locations of the vortices and domain walls and r(s)] have 
been eliminated, either by integrating them out or by 
setting them to their stationary values. For a numerical 
implementation of the latter procedure applied to super- 
conductivity in Sr2Ru04, see Ref. [15]. In the present 
paper, our aim is to proceed analytically, a task that is 
eased by our working in a particular limiting regime, an 
elaboration of the standard London limit that we term 
the "extended London limit" We remind the reader that 
the standard London limit amounts to assuming that the 
Ginzburg-Landau parameter k tends to infinity, which 
enforces the condition |r/| = 1. To pass to the extended 
London limit we make the additional assumption that 
the width of domain walls, which is controlled by the pa- 
rameter L, tends to zero. In this limit, the domain wall 
becomes vanishingly thin, compared with the penetration 
depth. 

We begin with the Ginzburg-Landau free energy, 
Eq. Q, and first pass to the standard London limit. 
From Eq. we see that in this limit the order- 

parameter amplitude \r]\ is energetically prohibited from 
departing from unity; inserting the corresponding form 
of the order parameter [i.e., Eq. ( [5a] ) but with \r]\ — 1] 
into Eq. Q, making the definition {AQ)ai '■— (Va^ — 
^a, Va7, Va/3)i, and dropping constant terms arising 
from the potential terms, we arrive at the unconven- 
tional superconducting free energy Fi^ corresponding to 
the London free energy for conventional superconductiv- 
ity, i.e.. 



F^[e,^,(3,A,H] 



d r<j -Kabcd^bidj (A6)ai (Ae)cj 



+ ^cos2/3+i|(Vx A)-Jf|2|, (17) 

^ai := Rli,iiIb^,-E^^,iM^J,{cosP,isml3)c, 



The aim of this section is to begin with the description 
of the superconducting system in terms of a Ginzburg- 
Landau free-energy functional dependent on the order- 
parameter field, and to derive from it a reduced descrip- 
tion in terms of the vortex and domain-wall densities 
and r(s) along domain walls. In this reduced descrip- 



where repeated indices i, j, . . . are summed from 1 to 3. 
Because in this free energy, the coefficient S is contracted 
with a tensor that is symmetric under time reversal, we 
have adopted a form for 5 that is manifestly symmetric 
under time reversal. In general, the supercurrent density 
J(r) is given by —SFsc/5A{r) and, continuing within the 
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London limit, we see that it has the form 

Ja-= - sin/3 Vq7 
1 



(18a) 



+ -(2/icos/3K6 + TXL)V6/3, (18b) 



1 



gab-= lab + l^TCQSPM^^ 



7— 7r/4 



(18d) 



Note the occurrence of the unconventional contribution 
JT" to the supercurrent, which includes currents that are 
localized near domain walls |19| . This contribution is 
manifestly odd under time reversal (which is evident be- 
cause each term is odd in /?). 

We now proceed to take the extended London limit, in 
which domain walls are controlled to be thin compared 
with the penetration depth. We begin by noting that 
the term arising from /p that remains in the free-energy 
density in the London limit is cos^ (/3) / 8L^ , and that this 
term contributes positively for any value of of /3 ^ ±7r/2. 
In particular, for a domain wall, across which /3 varies 
from 7r/2 to — 7r/2, the balancing, in equilibrium, of this 
potential term against contributions to the free energy 
that result from gradients in /3 would produce a spatial 
configuration in which /? changes from 7r/2 to — 7r/2 over a 
lengthscale (i.e., the domain- wall width) proportional to 
L. Thus, in the limit L — >■ 0, the widths of domain walls 
are controlled to be arbitrarily small, compared with the 
penetration depth (which, we remind the reader, we have 
chosen to set the unit for lengths). 

This extension of the London limit results in useful 
simplifications. First, as the domain walls are arbitrarily 
thin, regions in which (3 is uniform and equal to ±7r/2 
dominate, areally. Thus, terms proportional to cos/3 or 
sin/3 become or sgn/3 respectively. [Note that sgn/3 is 
the unit step function, taking the values 1 (or —1) for 
regions of positive (negative) chirality i.e., /3 > (or /3 < 
0)]. As an explicit example, the term in the superfluid 
density tensor g proportional to cos j3 can be neglected in 
the extended London limit, and thus we may make the 
replacement g ^ I. Physically, this means that, even 
in the presence of domain walls, the in-plane Meissner 
response is isotropic. 

A second useful simplification that arises in the ex- 
tended London limit is that it enables us to express con- 
tributions to the free energy and supercurrent involving 
gradients of /3 in terms the domain- wall density p^^ , de- 
fined in Eq. (13 1. Using Eq. (18) we can thus, e.g., write 



J=V0-A + J, 

J:= -sgn/3 V7 + {fiE + ^tM^^" 



(19a) 
(19b) 



For the sake of compactness, here and elsewhere we use 
the notation A^'"'-*-' • p^^ as shorthand for J dsA4^^''^ ■ 
n{s)S{r- R{s)). 



Having discussed how, in the extended London limit, 
the spatial variation of /3 is fully incorporated via the lo- 
cations of the domain walls {i?„(s)}, we continue with 
our goal of constructing an effective free energy by elimi- 
nating all degrees of freedom associated with excitations 
of the order parameter around the state of minimum free 
energy within a fixed topological sector. With this in 
mind, our next step is to eliminate the non-topological 
variations in the 6 field. 

As the superconducting order may possess vortices, 9 
is not, in general, a single- valued function of position, 
and therefore it may exhibit singular behavior (i.e., at the 
cores of vortices). Our initial strategy for eliminating the 
non-topological variations of 9 is to decompose it into two 
components: 9 — ^sm + ^vj where 9^^ is a smooth, single- 
valued part, and 9v is the part that accounts for any 
vortex singularities. This separation is not unique, but 
we shall see, after eliminating 9sin from the free energy 
by setting it equal to its stationary value 0sm, that the 
resulting free energy is — for any fixed choice of topolog- 
ical variables, such as vortex positions and strengths — 
independent of any particular choice of decomposition. 
To implement this elimination of 0sm we need only con- 
sider the terms in the free-energy density associated with 
the kinetic energy of the supercurrents (i.e., associated 
with J^), as other terms do not depend on 9; in the ex- 
tended London limit the free energy Fj constructed from 
these terms is given by 



Fj^ I d\;:\V9, 



-A + J\ 



(20) 



Stationarity of this expression with respect to 9sm reads 



= V ■iV9,-A + J), 



(21) 



and, by using the Green function for the Laplace oper- 
ator in two dimensions [i.e., G(r) = — ^ln|r|, obeying 
— V^G'(r) — 6{r)], the stationary solution 9sni can be ex- 
pressed as 

9sm{r')^ Jd\Gir'^r)V-{V9,-A + j){r). (22) 



By inserting 0sm into Eq. ( 20 ) , we arrive at the following 
form for the free energy: 

Fj^lJ d\\j^\\ (23a) 

Jj(r) := I d\' 5{v - r') - VaG(r - r') V',) 

X (Vfc^v - Afc + Jb) (23b) 

- j d\' EabVb G(r - r') E.^K (v^ 9^ - + Jo) , 

where we have used the elementary result EabEcd — 
lachd — ladhc and the defining equation obeyed by G. 
The procedure of minimizing Fj with respect to 0sm can 
be described, physically, as compensating for any source 
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of longitudinal currents (i.e., current- flows that build up 
at some location) or, equivalently, as a projection on to 
the subspace of transverse currents. This construction 
brings to the fore the vorticial content of the transverse 
supercurrent, which arises both from vortices and domain 
walls. Specifically, one can identify the vorticity W via 



Under the constraint that it be evaluated at the station- 
ary value of A, this form for Fj has the same value as the 
one given in Eq. ([23]). We note, in passing, that the dual 
free energy Fj [A] depends explicitly on the local value of 
the perpendicular magnetic field B {— Eab^ a^b)- Thus, 
the full expression for the free energy in the extended 
London limit, which also includes the magnetic field en- 



W = EabVa{Jb + Ab) = EabVa{VbO^ + Jb) ■ (24) ergy / (Prl{B-Hf, depends on B locally. This locality 



Owing to the unconventional contribution to the super- 
current iT", the vorticity W in unconventional supercon- 
ductivity in the extended London limit comprises both a 
vortex term Wv, which is common also to conventional 
superconductivity and is proportional to the total vor- 
tex density pv, and a domain- wall term Wdw, which is 
proportional to the domain- wall density pdw'- 



1 



(25a) 



:=27r(^-(l + sgn/3)p+(r) 

+ i(l-sgn/3)p7(r) 
Wdw := (/(s) n(s) + d{s) ■ V) • pdw, (25c) 



(25b) 



where 



-2 9,r(s), 



d{s) :— —fil 



. jj2(0(s)-r(s))^ 



(26a) 
(26b) 



Several points are worth noting here. First, is a 
weighted sum of the vortex densities in the chiral do- 
mains, which makes evident the fact that only those sin- 
gularities of {9~) that are located in the positive- 
chirality (negative-chirality) domain are associated with 
local vorticity. Second, via Eq. (25c), we see that the do- 



main wall vorticity Wd can be expressed as a sum of two 
contributions: a "monopole" contribution of strength /, 
which determines the net magnetic flux penetrating the 
superconductor; and a "dipole" contribution of strength 
d, which is generated by currents that flow along domain 
wall cores but do not create net flux through the super- 
conductor. Third, within this extended London limit, 
the monopole and dipole contributions are expressible in 
terms of the topological degrees of freedom r(s) and (t>{s). 
We remind the reader that (/()(s) is determined from the 
trajectory of a domain wall R{s) via Eq. (12). 



The final step in deriving the reduced free energy is to 
eliminate the vector potential A. Although it is possible 
to proceed directly, using Eq. (23 1 (see Appendix [b]), the 



fact that the current J in Eq. (23) is determined via 



a nonlocal expression makes it more efficient to apply 
an alternative, 'dual' approach, which uses a Hubbard- 
Stratonovich transformation of the nonlocal kernel via an 
auxiliary field A; see, e.g., Ref. [321 [SO]- The resulting, 
dual expression for the free energy Fj is thus given by 



Fj[A]^ J d\ l^~^\VA\^+A{W-B) 



(27) 



renders simple the task of identifying the stationary value 
of B. Eliminating B by setting it equal to its stationary 
value we arrive at the following form for the extended 
London limit of the free energy: 



i^EL=y"rfV |-iA(-V2 + l)A+A(T4^-H)+/corc| 



(28) 



In this form, the first two terms, which together account 
for the kinetic energy of the supercurrent and the mag- 
netic field energy, have the virtue of being local and de- 
termined via W (i.e., the vorticity of the supercurrent). 
The remaining contributions to the free energy given by 
Eq. (17 1 are accounted for via /core, which is associated 



with the core energy of the domain walls and are negli- 
gibly small in regions lying beyond a distance of a few 
wall widths L from a domain wall. An explicit expres- 
sion for /core in terms of the fields 7 and /3 is given in 
Appendix [A] [M] . Thus, in the neighborhood of the ex- 
tended London limit, in which L becomes small (but re- 
mains non-zero), the domain- wall energy Jd^r /core can 
be expressed in terms of an energy per unit domain-wall 
length Ecore, which depends locally upon on T{s) (i.e., 7 
evaluated on the domain wall) together with the shape 
of the domain wall [(e.g., via 0(s))], along with their ar- 
clength derivatives: 



dV/c, 



E 



ds Ecoic (r„(s), . . . ; (f>n{s), ...). 



(29) 

We are now in the position to complete our derivation 
of the reduced free energy i^EL in the extended London 
limit, reduced in the sense that it depends only on the 
external applied magnetic field and the topological vari- 
ables via the vorticity W and and domain-wall core en- 
ergy density /core- Upon eliminating A from Eq. (28 1, 
Fel becomes 



EL 



rf'r/c, 



(30) 



An 



{W{r)~H{r))Ko{\r-r'\){Wir')-H{r')), 



where Kq is a modified Bessel function of the second 
kind. A virtue of the formulation that we have employed 
is that it enables the efficient calculation of the magnetic 
response of the superconductor in the extended London 
limit, via the thermodynamic relation 



AfEL(r) 



(5F, 



EL 



SH{r) 



Ko{\r-r'\){W(r')-H{r')). 

(31) 
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One can also use the Ampere-Maxwell law to determine 
the spatial distribution of equilibrium supercurrents in 
the this limit, which gives Jel = E ■ VMel- 

As an initial illustration of this approach, we consider 
a straight domain wall, lying along the y-axis in infi- 
nite, three-dimensional superconductor. We assume that 
there is no applied magnetic field, i.e., H = 0. We fur- 
ther assume that the superconducting state is of positive 
(negative) chirality for a; < (x > 0), so that by the 
convention defined by Eq. ( [l2| we have R{s) = sy. As 
we show in the Appendix [Xj a variational analysis, based 
on an assumed form for the behavior of /3 transverse to a 
translationally invariant domain wall, suggests that the 
equilibrium value of F is (p. Assuming this to be case, we 
then find, from Eq. (26), that the domain- wall vorticity 
has no monopole part [i.e., /(s) = 0] but does have a 
dipole part, which is given by d = {^{■kt/A) — h)I [52] . 



Then, from Eq. (31 ), we find that magnetization and cur- 



rent densities vary with the spatial distance x from the 
domain wall as follows: 

Me^{x) = -i((Ai- (7rr/4))sgn(a;)e-l^l, (32a) 

JEL(a;) (ttt/A)) {6{x) - e-l-l/2) y. (32b) 

As, for this domain- wall configuration, the monopole con- 
tribution /(s) to the domain-wall vorticity is zero, the 
net magnetic flux (per unit length of domain wall) [e.g. 
the magnetic flux (per unit length of domain wall) inte- 
grated transversally] vanishes. The jump discontinuity 
in M{x) at a; = results from a supercurrent that flows 
along the domain- wall core. For the case of Sr2Ru04, 
we can use Eq. (32b I to estimate the magnitude of this 



current. In SI units the dimensionful current density J' 
is given in terms of its dimensionless counterpart J via 
J' = {2TTXfo/^o)J. Then, using J to compute the cur- 
rent passing through a narrow window bracketing the do- 
main wall, we arrive at the following expression for the 
dimensionful domain-wall current /' per Ru-0 layer: 

/' = '^^^S^Az f dx in - (^t/4)) 6ix) (33a) 



$0 



Az (/I - (ttt/A)) , 



(33b) 



where Az is the thickness of an Ru-0 layer. To find 
the numerical value of this current in Amps, we note 
that in SI units /o can be expressed in terms of the 
thermodynamic critical field He as 2k,^^iiqH'^. Using 
the parameter values appropriate for Sr2Ru04 [T], i.e., 
/ifliJc = 0.023 T, A(0) = 0.15 /im, n = 2.3, and Az = 
1.2 nm, we arrive at the following estimate for the cur- 
rent: (/i - (7rr/4)) X 1.3 x IQ-^ A per Ru-0 layer flow- 
ing along a domain-wall core. For this result to match 
previously made theoretical estimates (see Ref. [I5ll53]). 
one would need to have the material perameters obey 
(m-(W4)) 

In the next section, we extend our discussion to cope 
with situations lying beyond straight domain walls, thus 



allowing the domain walls to have bends. As part of 
this discussion, we employ the reduced description of the 
superconductor in the extended London limit derived in 
the present section to show that: (i) a net magnetic flux 
penetrates the superconductor near bends; and (ii) this 
flux is generically a nonintegral multiple of the supercon- 
ducting fiux quantum $o- 



V. MAGNETIC FLUX IN THE VICINITY OF A 
BEND IN A DOMAIN WALL 

In this section we derive the central result of this work, 
viz., that a bend in a domain wall is accompanied by 
a nonintegral amount of magnetic flux that penetrates 
the superconductor near the bend; the amount — which 
we term the bend flux — depends on the geometry of 
the bend. In the limit in which the in-plane crystalline 
anisotropy is negligible (i.e., the isotropic limit), the bend 
fiux is proportional to the angle through which the do- 
main wall bends. 

We derive the bend flux via two approaches. In the 
first, we analyze a bending domain wall via the effective 
theory of the topological variables, developed in Sec. |IV| 
We then consider an alternative derivation, which, in the 
isotropic limit, yields the bend flux quite generally, with- 
out reliance on the assumption of either the standard or 
the extended London limit, or even on the validity of the 
Ginzburg-Landau expansion of the free energy. We end 
this section by considering modifications of the isotropic- 
limit bend flux result that would arise in settings of other 
pairing symmetries and/or tetragonal or hexagonal de- 
partures from the limit of crystalline isotropy. 



A. Comparison with an a spatially extended 
Josephson junction 

Before establishing the existence of bend fiux, we give 
a discussion of the the essential differences between, on 
the one hand, a system comprising a domain wall and the 
superconducting regions of opposing chirality separated 
by it, and, on the other hand, a system of a spatially 
extended Josephson junction and two regions of conven- 
tional superconductivity coupled by it. For the extended 
Josephson-junction system it is possible to define a vari- 
able analogous to the domain wall variable r(s), i.e., the 
local value rj(s) := (-6li(s) + 6i2(s))/2 of (haff of the) the 
difference between the phases 9i{s) and ^2(5) of the su- 
perconducting regions that lie on either side of the junc- 
tion. The important distinction between F and Fj is that 
whereas Fj transforms trivially under in-plane rotations, 
F transforms nontrivially. 

This observation has important implications, if we 
compare the local energy of a domain wall i?dw with the 
Josephson energy Ej of a extended Josephson junction. 
In particular, for the extended Josephson junction, Ej is 
a periodic function of Fj alone. For the domain wall sys- 
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FIG. 4. (a) "Dipole" and (b) "monopole" contributions to the magnetic field associated with a bending domain wall (oriented 
black line) for a bend angle of 7r/3, considered in the extended London limit [see Eqs. ( 26 1 and (31 1]. The 2-axis magnetic field 
M(r) is plotted as a function of position (color scale and contour lines). For this example, the Ginzburg-Landau parameters 
are taken to have the values = 0.2 and t — 1. The penetration depth defines the unit lengthscale. The chirality is positive to 
the left of the domain wall and negative to the right of it, so that, via Eq. ( 12 1, we see that the domain wall has the indicated 
orientation. Only the monopole contribution produces a net magnetic flux through the system, (color online) 



tern, on the other hand, in order to preserve the 80(2)^ 
invariance of the free energy, the local energy-density of 
the domain wall must be a periodic function of the dif- 
ference T — (j) [in which (jj is determined by the local di- 
rection of the domain wall normal; see Eq. (12|]. Thus, 
because they have distinct values of 0, two segments of 
straight domain wall separated by a bend will generically 
have distinct equilibrium values of F. This stands in con- 
trast with the case of the spatially extended Josephson 
junction with a bend, the equilibrium value of Fj being 
independent of position along the junction. As we shall 
now see, the bend flux originates in this variation of the 
equilibrium value of F on either side of a bend. 



B. Bend flux in terms of topological variables 

We now turn to the derivation of the bend flux within 
the special context of the effective theory for topologi- 
cal variables, developed in the previous section. Part of 
the utility of this effective theory is that it allows for 
an efficient calculation of the magnetic response of the 
superconductor, given a conflguration of the topological 
variables, viz., the position and strength of the vortices, 
the positions of the domain-wall lines, and the value of 
F(s) along each such line. Thus, our approach will be 
to consider a specified configuration of topological vari- 
ables without vortices but with a single, fixed domain 



wall having a bend and a specified form for F(s) along it, 
and then to employ Eq. (31) in order to determine the 



corresponding magnetization density. 

We define the position of the domain wall using three 
line-segments: an arc of radians and unit radius of cur- 
vature, and two straight segments that continue tangen- 
tially from each of the end-points of the arc (see Fig. |4| . 
Given this particular geometry, we say that the result- 
ing domain wall has a bend angle of in it. Our next 
assumption concerns the form of T{s). In Appendix [A| 
we give a variational analysis that suggests that, for a 
straight domain wall, the equilibrium value of F is 0. 
To generalize to the situation in which the domain wall 
bends, we assume that F follows the local direction of the 
domain wall 'adiabaticallyj i.e., F(s) = 4>{s). In this case, 
because dsT(s) is not everywhere zero a monopole contri- 



bution to the domain- wall vorticity arises [see Eqs. (25c 
and (26a|], in addition to the dipole contribution. Fig- 



ure |4 shows both the monopole and dipole contributions 
to the magnetic field, evaluated using Eq. (31 1. 

Next, we determine the total flux $var,bcnd associated 
with the bent domain wall furnished by this variational 
calculation. To do this, we integrate the total magnetic 
field through a large circular disc centered at the vertex 
formed by the extrapolation of the straight-line segments, 
so that the straight-line segments lie radial to the disc. 
In the limit that the disc radius is much larger than the 
penetration depth, we find that the dipole contribution 
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ai- 
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^ 



f [ f 5. A bent domain wall (oriented black line) separat- 
ing two domains of opposite chirality. En route to deriving 
Eq. ( 40 1 , which expresses the the flux $ through the area 
bounded by the contour C in terms of the bend angle Q, the 
circulation of the unconventional superfluid velocity V [see 
Eq.( |35| l] around C is shown to be zero. The contour C is as- 
sumed to be many penetration depths away from the region 
where the domain wall bends, (color online) 



to *&var,bend tends to zero, whereas the monopole contri- 
bution is nonzero, tending to the following total flux: 



(frM{x) = 2 / dT{s) = 26, (34) 



var,bcnd 



i.e., the net flux is simply given by twice the bend angle, 
regardless of how r(s) interpolates between its limiting 
values far from the bend. In particular, for case shown 
in Fig. |4] (i.e., for 8 = 7''/3) the bend flux is 27r/3, i.e., 
the dimensionful value is $o/3, which is a nonintegral 
multiple of the flux quantum. 



C. General analysis for the bend flux 

In the remainder of this section we consider a more 
general context in which the existence of bend flux can 
be demonstrated. In particular, we need not employ the 
extended London limit, nor assume that the supercon- 
ductor is in the Ginzburg-Landau regime. Rather, the 
central assumption is that the superconducting order has 
the following essential feature: in regions in which the 
chirality is maximal, local S0{2)z rotations of the super- 
conducting order can equally well be accomplished via 
U(l) gauge transformations, so that the local transfor- 
mation (a, — sgn(/3) a) e U(l) x S0(2)z acts trivially on 
the superconducting order parameter. If this assumption 
holds then, provided the amplitude of the superconduct- 
ing order is spatially homogeneous away from the domain 
wall, the unconventional superfluid velocity V, defined 



V := V6I - sgn(/?) V7 - A, 



(35) 



tends to zero within a maximally chiral region, as a result 
of the Meissner effect. 

In deriving the bend flux we use the following con- 
struction to define the geometry of the domain wall. We 
consider a single domain wall that is fixed to pass through 
three points: the origin O, and two other points P and 
Q; then, we fix the angle ZPOQ — tt + Q where O is 
bend angle of the domain wall. We take the orientation 
of the domain wall, as defined by Eq. (12), to run from 
P to Q, and we let the points P and Q tend to infin- 
ity. The goal of the calculation is then to determine the 
net magnetic flux penetrating the superconductor in the 
vicinity of the domain wall bend. The key quantity that 
we use is the circulation of V around a closed contour 
C encircling the domain wall bend at a distance of many 
penetration depths (see Fig. [5]). Care is needed in select- 
ing the contour C because, in equilibrium, even for zero 
applied magnetic field, a domain wall may not be trans- 
lationally invariant |54j . However, as the underlying free 
energy is local [as is manifest in Eq. (28)] and transla- 
tionally and rotationally invariant; it is always possible 
to choose two geometrically congruent regions, A and B, 
each straddling the domain wall but located on opposite 
sides of the domain-wall bend, in which the equilibrium 
spatial configurations of the superconducting order in the 
regions (A and B) are related to one another via a rigid 
rotation and translation (see Fig. 5]). Once a pair of such 
regions has been identified, we choose the contour C to 
cross the regions (and hence the domain wall) on locally 
identical paths 61 — ?> 63 and ai — >■ 03 (i.e., on paths that 
are related by the same rotation and translation as the 
regions are). As a result, the following equality between 
line integrals holds: 



I dr-V^ j 



dr-V. 



(36) 



By using this result, and observing that V = away 
from the domain wall, we see that the circulation of Y 
around the closed contour C is zero, i.e.. 



dv-V ^Q. 



(37) 



The next step in the derivation is to consider the con- 
tour (C~), which begins at the point a2 (^2) and fol- 
lows C through the positive- (negative-) chirality domain 
to the point 62 (02)- The line-integrals of (V -I- A) along 
C"*" and C~ respectively measure the change in the phase 
of the order parameter in the positive (negative) region 
from 02 to 62 (and from &2 to 02). Thus, ag ain using 
the linear relation 7 = {-9+ + 9^)/2 [i.e., Eq. ([lO|], we 
see that the change in F from the point 62 to the point 
02 [i.e., AF := F(a2) — F(62)], is given by the following 
formula: 



AF = 



\l d9- =^-<f dr-iV + A). (38) 
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We now examine in more detail the equilibrium value of 
Ar. As discussed in Sec. |V A as a consequence of the 
rotational invariance of the underlying free energy, the 
energy (per unit arclength) _Edw of the domain wall must 
be a periodic function of the combination T{s) — (t>{s), in 
which (/)(s) continues to characterize the local direction 
normal to the domain wall. Furthermore, as — up to a 
global phase — the configurations having T and F + tt are 
equivalent, the dependence of i?dw on F(s) — 4>{s) has 
period tt. Importantly, we make the following additional 
assumption, viz., that the dependence of i?dw on F(s) — 
(j){s) has a single minimum per period. 

We now observe that, by construction, region A is 
rotated by an angle Q relative to B (using the con- 
vention that positive rotations are measured counter- 
clockwise, relative to the domain-wall orientation) as a 
result A(j) :— (t>{ai) — 4>[a2) — 6. Thus, with these as- 
sumptions the equilibrium value of AF is equal to the 
bend angle Q, up to an integer multiple of tt, i.e.. 



AF = e 



(39) 



Combining Eqs. ( 37|38 39), and defining $bGnd to be the 
bend flirx (i.e., net flux through the surface defined by 
the contour C) we arrive at the result that 



bend 



{{e/TT)+n) $0- 



(40) 



Because the bend flux, in the rotationally invariant limit, 
can evidently be an arbitrary fraction of the flux quan- 
tum, this result is a manifestation of the general result 
that broken time-reversal invariance allows for nonquan- 
tized amounts of flux to penetrate a superconductor, as 
predicted on general grounds in Refs. [2^ [2^ fSBHSUlISS] . 
Moreover, because F(s) need not stay locked, relative to 
the local domain wall orientation (e.g, at the bend), or 
owing to the presence of vortices in either or both of the 
chiral domains, it makes sense that $bcnd be determined 
only modulo $o f56l. 



D. Bend fluxes for other pairing and crystalline 
symmetries 

In this paper, we have assumed that the superconduct- 
ing order transforms as one particular representation of 
80(2)^. We now obtain the generalization of the for- 
mula for the bend flux, Eq. (40), that remains valid for 



arbitrary irreducible representations, which can be in- 
dexed in terms of an integer m (see, e.g., Ref. [S7]). For 
brevity's sake, we refer to the m = I case as p-wave 
(which is the case focused on in this paper), and the 
m = 2 case as d-wave. En route to generalizing Eq. (40) 



to arbitrary m, we assume that transformations of the 
form (ma, — sgn(/3) a) G U(l) x 80(2)^ act trivially on a 
uniform, maximally chiral phase. Under this assumption, 
and repeating the line of argument given in 8ec. |V C] mu- 
tatis mutandis, the bend flux formula becomes 



Another version of Eq. (40 1 results when we address 
the setting of tetragonal D4^h symmetry (which is, of 
course, discrete). In this case, the argument given in 



{{niQ/n) + n) $0. 



8ec. |Vq leading to Eq. (|40| holds only for 9 = ±tt/2, 
for which the minimum net flux through the domain wall 
bend in the p-wave case is $o/2; this is distinct from a 
conventional vortex, for which the net flux is $0 . In con- 
trast, for the d-wave case and 8 = ±7r/2, a net flux of 
n$o (with n integral) penetrates the bend. As a last ob- 
servation, we note that for p-wave pairing and Dq^ sym- 
metry and a domain- wall bend angle of 7r/3, the smallest 
positive net flux accompanying the bend is <I'o/3, whereas 
the smallest negative net flux accompanying it is — 2'I>o/3. 

Now that we have established that in, various settings, 
one anticipates that a bent domain wall is accompanied 
by nonintegral net magnetic flux, we shall, in the next 
section, explore a range of experimental situations in 
which such effects might be observable. 



VI. EXPERIMENTAL IMPLICATIONS 

We now describe three experimental scenarios in which 
it may prove possible to observe, via scanning magnetic 
microscopy, the phenomenon of sample penetration by 
nonintegral net magnetic flux associated with bends in 
walls separating regions of opposing chiral superconduct- 
ing order. These scenarios are depicted schematically in 
Figs.[6j[7j and[8j Augmenting the bend flux phenomenon 
discussed in the present work, it is known that domain 
walls are expected to produce Ampere magnetic flelds, re- 
sulting from currents that flow along the core of a domain 
walls; see Refs. [T21[TS]. We emphasize that, even if there 
were a specific, microscopic reason for the magnitude of 
such currents to be reduced (cf., e.g., Refs. [T H \58 \ [59 ] ). 
e.g., below currently detectable levels [TS], such a reduc- 
tion would not affect the existence or magnitude of the 
bend fluxes discussed here. Thus, it is perhaps useful to 
regard bend fluxes as providing robust magnetic signa- 
tures of domain walls, and hence the form of supercon- 
ductivity that spontaneously lacks time-reversal symme- 
try. 

As a first example, consider a domain wall that inter- 
sects a physical surface of the superconducting system, 
the surface being oriented perpendicular to the z axis. 
In this paper we are neglecting effects resulting from the 
finite height above the sample surface at which magnetic 
fields would typically be detected. (For a discussion of 
such effects see, e.g., Ref. [5^.) In addition, we envision 
domain walls to be pinned at generically located sites, 
e.g., by impurities. In the limit in which the bulk terms 
in the free energy that break 80(2)^ symmetry are small 
(as can hold occur near T^.), the spatial arrangement of 
these pinning sites predominates in determining the bend 
angles that characterize a domain wall as it traverses the 
sample. Assuming that these pinning sites are spaced 
further apart than the penetration depth, Eq. ( 40 ) indi- 



(41) cates that these pinning locations would show up in scan- 
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$2 = -262 




FIG. 6. Schematic depiction of a chiral domain wall running 
along an ab-face of a superconductor and pinned to various 
locations. The wall is indicated by the oriented black line 
[the orientation is defined by Eq. ( 12 l]. We assume that devi- 
ations from S0(2)z symmetry are sufficiently small that the 
arrangement of the pinning sites determines the path of the 
domain wall. Orange dots denote pinning sites. Near them, 
the domain wall bends and flux penetrates the superconduc- 
tor. The specified positive angles {Si} express the geometry 
of the bends, as indicated. The bend fluxes are then deter- 
mined via E q. (|40[ ). For each bend, we have chosen the value 
of n in Eq. ( |40[ ) to give the corresponding bend fluxes {"I>i} 
the smallest possible magnitudes. The orientation of the flux 
accompanying each bend is indicated via a dot (up) or a cross 
(down), (color online) 



FIG. 7. Schematic depiction of an array of domain walls (ori- 
ented black lines), each intersecting an ab-face of a supercon- 
ductor. D4h deviations from SO(2)z symmetry are assumed 
to be large enough to pin domain walls to lie along specific 
crystallographic directions, and each domain wall is assumed 
to have a 7r/2 bend. Integral- flux vortices (intermediate size 
orange dots) penetrate the superconductor along the straight 
sections of the domain walls. At each bend, a bend flux pen- 
etrates the superconductor (large and small orange dots) and 
is fixed, via Eq. (401, to be a half-integer multiple of the flux 
quantum. All regions of localized flux are shown as if they 
had the same sign, as would be energetically favorable in the 
presence of a magnetic field applied along the z axis, (color 
online) 



ning magnetometery as local regions of nonquantized flux 
penetrating the superconductor (see Fig. [6]) . 

We now outline a scenario specific to Sr2Ru04. In 
both zero and nonzero in-plane magnetic fields, scan- 
ning magnetic imaging of Sr2Ru04 shows that vortices 
arrange themselves in line- like structures [161 IHl EB l£2] ■ 
One of the possible scenarios put forth to explain these 
structures is that the line-like structures are due to the 
binding of vortices to a parallel array of chiral domain 
walls 26J. However, to date, the line-like structures have 
not exhibited characteristics that would uniquely identify 
them as domain walls because, to within experimental 
uncertainty, the vortices (i.e., the local regions of pene- 
trating magnetic field) were observed to have total fluxes 
that were integer multiples of $oj and Ampere magnetic 
fields along the line-like features were not observed. The 
phenomenon of bend flux provides an additional route for 
determining whether the observed line-like structures are 
indeed associated with domain walls. If it proves possible 
to prepare a sample (e.g., via a fleld-sweep procedure) so 
that the line-like features are bent then, if the line-like 
structures do indeed correspond to domain walls, bends 
would be accompanied by a nonintegral flux penetrating 
the superconductor (see Fig. [t]) . The observation of non- 
integral bend flux at a 7r/2 bend would provide further 
confirmation of p-wave pairing in Sr2Ru04 because, as 
noted in Sec. |VD[ d-wave pairing would produce integer 
bend flux. However, the fact that Refs. [151 EH El] do 



not report regions of nonintegral localized flux suggests 
that, in these experiments, if there are domain walls then 
they are aligned in parallel arrays, and thus are not bent. 

A further consequence of domain walls should be ev- 
ident in annular rings of broken time-reversal symme- 
try superconductors. As the analysis leading to Eq. (40) 



is local only to the contour C, and does not require in- 
spection of the superconductivity near the domain-wall 
bend itself, it can be generalized to the case in which the 
bend is replaced by a hole in the superconductivity (see 
Fig. |8]) . Recently, evidence for half- integer fluxoid behav- 
ior has been obtained in experiments on mesoscopic rings 
of Sr2Ru04 using cantilever torque magnetometry [S]. 
However, in those experiments the half-integer fluxoid 
behavior was found to be accompanied by a small, ro- 
tationally invariant, in-plane component of the magneti- 
zation, and we are not aware of any reason why such a 
magnetization component would arise in the context of 
domain walls. 



VII. CONCLUDING REMARKS 

We have analyzed the properties of unconventional su- 
perconductors in which the superconducting state spon- 
taneously breaks time-reversal symmetry and thus have 
the potential to exhibit domain walls that separate re- 
gions of opposing order-parameter chirality. By employ- 
ing an extension of the well-known London limit of the su- 
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$ = $o/2 



+ 



FIG. 8. Schematic depiction of an annular sample (shaded 
gray) the thickness and hieght of which is larger that the pen- 
etration depth. The annulus is crossed radially by a pair of 
domain walls oriented at 7r/2 relative to one another. Tetrag- 
onal D4h terms in the free energy are assumed to be large 
enough to pin the direction of the domain walls to crystallo- 
graphically to lie along the specified directions. The minimum 
magnitude of the total flux through the hole would then be 
"l?o/2. (color online) 



from its intrinsic interest, resolving this question would, 
inter alia, be valuable in assessing the utility of mate- 
rials such as these for exhibiting nonabelian phases and 
Majorana modes and, hence, the robustness with respect 
to decoherence that could prove useful for quantum in- 
formation processing purposes. 
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Appendix A: Free energy of a translationally 
invariant domain wall 



perconducting state, we have formulated an effective the- 
ory in terms of the topological variables that describe vor- 
tices and domain walls of the order parameter. We have 
used this effective-theory formulation to show that local- 
ized near any bend in a domain wall through an angle 0, 
there is an associated net magnetic flux ((O/tt) -t- n) $o 
(for some integer n) — provided the system can be taken 
to be rotationally invariant, crystallographically, about 
the z axis. We have also shown that this result for the 
flux near a domain-wall bend holds more generally. Nei- 
ther the London limit nor the regime of validity of the 
Ginzburg-Landau theory are required. Rather it is suffi- 
cient for the following condition to hold: within regions 
of maximal chirality the two transformations, 80(2)^ ro- 
tations and U(l) gauge transformations of the supercon- 
duction order parameter, are degenerate transformations, 
in the sense that they have equivalent impacts on the 
state of the superconductivity. 

We have addressed the issue of the relaxation of the 
assumption of crystallographic rotational invariance, and 
its replacement by discrete rotational invariance. In this 
regime we have found that the result for the bend flux 
continues to hold, but only for specific values of the bend 
angle, that are determined by the crystalline symmetry. 
We have also sketched three candidate settings in which 
the interplay between chiral-domain-wall geometry and 
magnetic flux discussed in this paper might be observ- 
able, e.g., in experiments using scanned probe magnetic 
imaging. The analysis that we have presented may be of 
use in determining the existence and distribution of do- 
main walls in various superconducting materials such as 
Sr2Ru04, and may thus be of use in resolving the ques- 
tion of whether superconductivity in Sr2Ru04 does in- 
deed spontaneously break time-reversal symmetry. Aside 



In this appendix we use a variational approach to de- 
rive an estimate for the free energy per unit length .Edw of 
a translationally invariant domain wall, starting from the 
free energy _Fel (in the extended London limit), given in 
Eq. (30 1. (Similar calculations can be found throughout 
the literature; see, e.g., Refs. [H HI HH [H [SH [Ml IM]-) 



As we discussed in Sec. |IV[ ^el contains two contribu- 
tions: one, /core, due to the core energy of a domain 
wall, which we estimate variationally; and the other, the 
London term, that describes both the kinetic energy of 
supercurrents and the magnetic field energy. In order to 
express /core compactly, we define (01,012) = and 
thus /core become 



/c( 

T 



1. 



1 

8i2 



COS /?, 



(Ala) 



1 /4cos^/3 







TCOS/3 -y,(^/4) / 4 

8 I 1-2^ 



/i sin 2(3 



Eab E, 



(Alb) 



where repeated indices {a, 6, are summed from 1 to 
2. The approach taken in this appendix is to evaluate in- 



dependently the two contributions to Eq. ( 30 1, expressing 



separately the variational estimate for the core energy per 
unit length i?corc and the London energy per unit length 
i?L, and then to add these contributions to determine 
Edw 

To derive the variational estimate, we make the follow- 
ing assumptions for the spatial dependence of the 7 and 
/? fields transverse to the domain wall: we take 7 to be 
constant and equal to T, and we take P{x) to be equal 
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to I3e{x) := 2tan~-'^ tanh(x/2£), where £ is a variational 
parameter specifying the width of the domain waU. [To 
motivate the form /3f , we note that for (/i, t) = (0, 0) and 

7 constant, the term /core reduces to g| V/3p + ^^^ cos^ /3, 
and this form has the property of being stationary at 
(3e{x), provided £ = L.] 

By using the variational assumptions for 7 and /3 we 
obtain the following expressions for Score, which depends 
upon £ and T as well as the angle <j) [which specifies the 
direction n — (cos(/),sin(/)) normal to the domain wall]: 



E,orc{r,q^,£) 

~ 1\4J^ ^ 4 ~ Y ~ 16 



By minimizing i?coro with respect to £, we see that, in the 
extended London limit (for which L tends to zero), the 
value of £ that makes E^^ stationary, would also tend 
to zero, provided the other energetic contribution, 
does not force the stationary value of £ away from this 
result. To see that indeed Ej, does not do this, using the 
same varriational assumptions for 7 and /? we examine 



-Bl expressed as power series in £ to 0{£'^): 

i?L(r, = 

+ I_cos4(r-0)+O(^)y (A3) 

Combining the two terms, £^core and i?L, we arrive at the 
following variational expression for the free energy per 
unit length of a translationally invariant domain wall: 

Ed^{T,<j),£) = 

+ y cos4(r-0) + O(^)j. (A4) 

By minimizing E^w with respect to £^ and recalling that 
in the extended London limit L is small, we find the sta- 
tionary value of ^ to be proportional to L, consistent with 
the assumption, just made, that £ is also small in the ex- 
tended London limit. Then, by replacing £ by its station- 
ary value one obtains a value for _Edw having the following 
properties, some of which we make use in Sec. |v] (i) it 
depends on T and (j) only through the combination F — 
and is n periodic in this quantity; (ii) it is independent of 
fi (to leading order in L); and (iii) when t < n/A the val- 
ues of F that minimizes min^ i?dw are (p + nn (for integer 
n). 



Appendix B: Free energy in terms of topological variables for the case of conventional superconductivity 

To motivate the derivation of the efl^ective free energy for the topological variables given in Sec. |IV|of the main text. 



resulting in Eq. ( 30 1 , we review in detail how it would proceed in the simpler setting of conventional superconductivity, 
and without employing the dual approach. For a conventional superconductor, the order parameter is the complex 
scalar field ip{r). We assume that, in the absence of a magnetic field, the system is translationally and rotationally 
invariant, and we consider magnetic fields that are oriented along the z direction and states of the superconductivity 
that are homogeneous in the z direction. In addition, we work with dependent and independent variables that have 
been rendered dimensionless via the rescalings given in Sec. [TTl 

With these assumptions we begin this derivation with the Ginzburg-Landau free energy per unit length of sample 

A, H]:^\j d\ {\{V - zA)V)|^ + y (IV^I^ - 1)2 + |(V X A) - H\'^ . (Bl) 

In the London limit, in which k — > 00, the potential terms of this free energy fix the magnitude of "0 to be unity. 
Then, -0 can be parametrized via a U(l) phase field 6'(r), so that ^(r) — expi0(r). Making this replacement in the 
free energy, we obtain the London form of the free energy, i.e., 

(fr{\Ve-A\^ + \{VxA)-H\^). (B2) 

The two terms in the London free energy can be regarded as frustrating one another, energetically, as they impose 
competing demands on the A field. The first term favors the transverse (i.e., divergence- free) part of A to be zero, the 
field can compensating for any longitudinal (i.e., curl-free) part; in contrast, the second term favors the transverse 
part of A to be nonzero. 
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For Type II superconductors at magnetic fields above the lower critical field, a partial resolution to this frustration 
comes from the introduction of vortices, which alter the structure of the 6 field: 9 becomes multi-valued, and is 
singular in the cores of the vortices. In particular, the expression V0 is not curl free and, correspondingly, has a 
transverse part. 

To derive the effective free energy in terms of the appropriate topological variables (in this case, the density of 
vortices) one now decomposes the 9 field into a smooth, single-valued part ^sm and a part 9^ that contains the vortex 
singularities, so that 9 — 6*31,1 + ^?v Next, one seeks to eliminate 9sn\ from the free energy by setting it to the value that 
makes the free energy stationary. As the only term in the free energy that depends on ^sm is the one corresponding 
to the kinetic energy of the supercurrents [i.e., the former term in Eq. (B2)], for the issue of stationarity one need 
only consider this term. Expanding it and integrating by parts, gives 

rf'^ Q^?sm(-V2)0sm - ^sm V • - A) + ^ I V^?v - . (B3) 

Then, using the Green function for the Laplacian in two dimensions, which obeys — V^G(r) — 5{y) and reads G{y) = 
— (27r)^^ In |r|, one finds that at stationarity 9^^ is given by 

^sm(r') = -y"d'rG(r'-r)V- (V0v(r)- A(r)). (B4) 



j2 i2 / 

a r a r 



By inserting 9sm into Eq. (B3) and using the defining equation for G(r) to express the last term of Eq. (B3) in terms 
of G, one obtains for the kinetic energy of the supercurrent 

(Va(V,^?v - Aa){r)) G{r - r'){v[{V,9, - Aa){r')) 

+ ({Va6. ~ Aa){r)) ( - V^G{r r')) ((V^^^ - ^a)(r'))) • (B5) 
The two terms in this equation have similar structure, and integration by parts allows them to be expressed as 

i J (fr (fr' [lablcd - ladhc) (Va^v - Aa) (r) (v6V,G(r - r')) {^'^9, - A^) (r')) . (B6) 
Next, by using the elementary tensor identity 

lab led — lad Ibc = Eac (B7) 

and integrating by parts, the suprcurrent kinetic energy becomes 



r' (^Eab (v„(Vb0v - Aa){r)) G(r - r') E.a (v^(V',0v - A„)(r'))) . (B8) 



This form shows that the elimination of the smooth part of 9 creates a long-ranged interaction for the curl of V^v — A- 
This free energy can readily be shown to be e quiv alent to Eq. ( 23 1 , and thus to describe the kinetic energy of the 
transverse part of the supercurrent. Equation (B8I features the curl of the gradient of the multi- valued function 9^, 
which is a combination that isolates the (5- function contributions from the singularities in the vortex cores, so that 

K&VaV60v = 27rpv, (B9) 

where Pv(r) := ~ -^i^) defines the vortex density in terms of the vortex locations {Ru} and vorticity {q^}- In 

particular, one sees that owing to the vortices the gradient of 9 can posses a transverse part, and this can partially 
relieve the frustration of A inherent in the London free energy. 

To proceed further with the derivation of the effective free energy in terms of vortex variables, one now considers 
the full London free energy, Eq. (B2), which, in terms of the total magnetic field B = Eab^aAb, reads 

i j d?r(fr' (27rpv - B){r) G(r - r') (27rpv - B){r') + ]^j(fr{B- Hf. (BIO) 

Note that we have omitted a constant contribution resulting from the suppression of the magnitude of the order 
parameter within the core of each vortex, as it is negligibly small, relative to the kinetic and field energies, in the 
London limit. 
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The next step is to eliminate the magnetic field from the free energy by setting it to its stationary value B which, 
from Eq. (BIO), one sees is 



^(r) = J (fr' (g + (5) \v- r') [h{y') + 2n J cfr" G{r' - r" 



)p{r": 



(Bll) 



where (G + (5) ^(r — r') is the inverse of the kernel G(r — r') + i5(r — r'). It is convenient to adopt a schematic notation 
in which one suppresses integral signs and dependences on spatial variables, in which case the result for B reads 



B ^ {G + 6)-\H + 2nGp). 



(B12) 



Replacing B hy B in Eq. (B10| then yields the following expression for the free energy: 
1 



((27rp) (G - G(G + Sr^G) (27rp) + H {S - {G + 5)-^) H - H{G + 6)-^G{2t:p) - {2ttp)G{G + Sy^H) .(B13) 

It is straightforward to see that each of the four integral kernels in this formula is the Green function for the Helmholtz 
operator in two dimensions, which obeys (— + l)G{r) = (5(r), and is given by Q{r) — (27r)^^Xo(|r|), where Ko{x) 
is a modified Bessel function. To exemplify this one can apply the following elementary manipulations to the kernel 
of the first term: 



G - G(G + 6)-^ G = G ((5 - ( - V2(G + (5)) ^ G ((-V^ + S)g - ^ G {- V^) G ^ g. 



(B14) 



By similarly simplifying the remaining kernels in Eq. ( |B13[ ) one completes the derivation of the effective free energy 
in terms of the vortex density pv and the applied field iJ, arriving at the result 



j- J d\£r' (2^pv - f) (r) i^o(k - r'l) (2^p^ - h) (r'), 



(B15) 



which is the analog for conventional superconductivity of the unconventional superconductivity formula Eq. ( 30 ) . 

I 



[1] A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 
(2003). 

[2] T. M. Rice and M. Sigrist, J. Phys. Cond. Mat. , L643 
(1995). 

[3] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. 
Mao, Y. Mori, and Y. Maeno, Nature 396, 658 (1998). 

[4] K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 
306, 1151 (2004). 

[5] J. Jang, D. Ferguson, V. Vakaryuk, R. Budakian, 
S. Chung, P. Goldbart, and Y. Maeno, (2010, unpub- 
lished). 

[6] N. B. Kopnin and M. M. Salomaa, Phys. Rev. B 44, 9667 
(1991). 

[7] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
[8] A. Y. Kitaev, Ann. Phys. (NY) 303, 2 (2003). 
[9] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 
S. D. Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
[10] Y. Yoshioka and K. Miyake, J. Phys. Soc. Jpn. 78, 074701 
(2009). 

[11] S. Raghu, A. Kapitulnik, and S. A. Kivelson, Phys. Rev. 

Lett. 105, 136401 (2010). 
[12] G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, 

J. Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Y. M. 

Z. Q. Mao, H. Nakamura, and M. Sigrist, Nature 394, 

558 (1998). 

[13] J. Xia, Y. Maeno, P. T. Beyersdorf, M. M. Fejer, and 



A. Kapitulnik, Phys. Rev. Lett. 97, 167002 (2006). 
[14] C. Kalhn and A. J. Berlinsky, J. Phys. Cond. Mat. 21, 
164210 (2009). 

[15] M. Matsumoto and M. Sigrist, J. Phys. Soc. Jpn. 68, 994 
(1999). 

[16] P. G. Bjornsson, Y. Maeno, M. E. Huber, and K. A. 

Moler, Phys. Rev. B 72, 012504 (2005). 
[17] J. R. Kirtley, C. Kallin, C. W. Hicks, E. A. Kim, Y. Liu, 

K. A. Moler, Y. Maeno, and K. D. Nelson, Phys. Rev. 

B 76, 014526 (2007). 
[18] C. W. Hicks, J. R. Kirtley, T. M. Lippman, N. C. 

Koshnick, M. E. Huber, Y. Maeno, W. M. Yuhasz, M. B. 

Maple, and K. A. Moler, Phys. Rev. B 81, 214501 (2010). 
[19] G. Volovik and L. Gor'kov, JETP 61, 843 (1985). 
[20] F. London and H. London, Proc. R. Soc. Lond. A 149, 

71 (1935). 

[21] V. Ginzburg and L. Landau, Zh. Eksp. Teor. Fiz. 20, 
1064 (1950). 

[22] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. 

Rev. 108, 1175 (1957). 
[23] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
[24] V. P. Mineev and K. V. Samokhin, Introduction to 

Unconventional Superconductivity (Gordon and Breach, 

Amsterdam, 1999). 
[25] M. Sigrist, T. M. Rice, and K. Ueda, Phys. Rev. Lett. 

63, 1727 (1989). 



19 



[26 
[27 
[28 
[29 
[30 
[31 
[32 
[33 

[34; 

[35 
[36 

[37; 

[38 

[39 
[40^ 
[41 

[42 

[43; 



Nat. Acad. Sci. USA 97, 2431 



M. Sigrist and D. F. Agterberg, Prog. Theor. Phys. 102, 
965 (1999). 

E. Babaev, A. Sudbo, and N. W. Ashcroft, Nature 431, 
666 (2004). 

G. E. Volovik and L. P. Gor'kov, ,]ETP Lett. 39, 674 
(1984). 

G. E. Volovik, Proc. 
(2000). 

M. Sigrist, D. B. Bailey, and R. B. Laughlin, Phys. Rev. 
Lett. 74, 3249 (1995). 

R. Heeb and D. F. Agterberg, Phys. Rev. B 59, 7076 

(1999) . 

Z. Q. Mao, Y. Maeno, S. NishiZaki, T. Akima, and 
T. Ishiguro, Phys. Rev. Lett. 84, 991 (2000). 
K. Deguchi, Z. Q. Mao, H. Yaguchi, and Y. Maeno, Phys. 
Rev. Lett. 92, 047002 (2004). 

S. Kittaka, H. Taniguchi, S. Yonezawa, H. Yaguchi, and 
Y. Maeno, Phys. Rev. B 81, 180510 (2010). 
We do not need to specify whether the representation is 
or , for which the basis functions transform respec- 
tively as {XZ, YZ} or {X, Y}. The results of the present 
work apply to both cases. 

J.-X. Zhu, C. S. Ting, J. L. Shen, and Z. D. Wang, Phys. 
Rev. B 56, 14093 (1997). 

C. Bergemann, S. R. Julian, A. P. Mackenzie, 
S. NishiZaki, and Y. Maeno, Phys. Rev. Lett. 84, 2662 

(2000) . 

M. W. Haverkort, I. S. Elfimov, L. H. Tjeng, G. A. 
Sawatzky, and A. Damascelli, Phys. Rev. Lett. 101, 
026406 (2008). 

N. D. Mermin, Rev. Mod. Phys. 51, 591 (1979). 
R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 
W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 
B 22, 2099 (1980). 

Y. Ran, Y. Zhang, and A. Vishwanath, Nat. Phys. 5, 
298 (2009). 

This choice of parametrization is similar to that used in 
Ref. [26], in which the aforementioned additional struc- 
ture of the order is parametrized by the scalar fields a 



[51] 



and X via rj oc R 



V2 



3(X + V4),^sin(x + 7^/4))• 



[44] We note that for weak coupling and specular reflection, 
the condition at the sample boundary requires /3 = 
(see, e.g., Ref. [23]), and thus, for a finite sample, all re- 
gions can be considered to be surrounded by domain-wall 
loops. However, the coupling between F and the normal 
direction of either the surface or the domain wall is gener- 
ically different, which can result in distinct equilibrium 
orientations of F for each of these cases (see, e.g. [15|). 

[45] Y. A. Izyumov and V. M. Laptev, Phase Transitions 20, 
95 (1990). 

[46] T. A. Tokuyasu, D. W. Hess, and J. A. Sauls, Phys. Rev. 

B 41, 8891 (1990). 
[47] J. A. Sauls and M. Eschrig, N. Jo. Phys. 11. 
[48] M. Ichioka, Y. Matsunaga, and K. Machida, Phys. Rev. 

B 71, 172510 (2005). 
[49] A. M. Polyakov, Gauge Fields and Strings (Harwood 

Academic Publishers, New York, 1987). 
[50] A. Zee, Quantum Field Theory in a Nutshell (Princeton 

University Press, Princeton, 2003). 



[52] 



[53] 
[54] 



[55] 
[56] 



Here and elsewhere in this paper, we take into account 
the core energies of domain walls but not the core energies 
of vortices. Our justification for doing this is that, in the 
standard London limit, the energy cost of a vortex core 
is negligibly small, compared with the kinetic energy of 
the supercurrents and magnetic fields, whereas the core 
energy of a domain wall is not. 

In the case that the equilibrium value of F — </> is 7r/2, the 
dipole part of the vorticity is given by d = ( — (7rr/4) — 
fijl. Thus, as is shown in the referring paragraph, the 
part of the domain-wall current proportional to r can 
fiow in either direction, for a given pattern of chirality, 
depending on the value of F — </!>. 

H.-J. Kwon, V. M. Yakovenko, and K. Sengupta, Syn- 
thetic Metals 133-134, 27 (2003). 

For example, as noted in Ref. ^5], in the limit of large 
domain-wall currents, vortices may be stabilized along 
the domain wall. Such vortices would then spoil the trans- 
lational invariance of the domain wall. Although we know 
of no experimental evidence for such an effect, there is 
evidence from Josephson-j unction tunneling experiments 
on Sr2Ru04 [65] for the absence of translational invari- 
ance (i.e., for the presence of multiple domains) along 
sample boundaries, which can be regarded in certain re- 
spects as analogous to domain walls. 
V. B. Geshkenbein, A. I. Larkin, and A. Barone, Phys. 
Rev. B 36, 235 (1987). 

As a particular case of Eq. (401, one can consider a 
straight domain wall. Then, Eqs. (39 1 and (40 1 imply 
that the topologically stable, localized solitons in F(s) 
along a domain wall would obey AF = tt, and that each 
is associated with a fiux $o- However, if the dependence 
of Sdw on F(s) — (j){s) should have multiple minima per 
TT period, then there could be topologically stable soli- 
tons in F(s) along a straight domain wall, each having 
AF 7^ TT and connected with nonquantized amounts of 
fiux (see Refs. [251 126;, 6 4 ). In this case, in addition to 
the bend flux of Eq. ( 40 1 , the flux associated with a bent 



[57] 



[58] 



[59] 

[60] 
[61] 



[62] 
[63] 
[64] 
[65] 



domain wall may have a further contribution. 
L. D. Landau and E. M. Lifshitz, Quantum Mechanics : 
Non-relativistic Theory, 3rd ed. (Pergamon Press, New 
York, 1977). 

A. J. Leggett, Quantum Liquids : Bose Condensation 
and Cooper Pairing in Condensed- Matter Systems (Ox- 
ford University Press, New York, 2006). 
P. E. C. Ashby and C. KaUin, Phys. Rev. B 79, 224509 
(2009). 

H. Bluhm, Phys. Rev. B 76, 144507 (2007). 
V. O. Dolocan, C. Veauvy, F. Servant, P. Lejay, K. Has- 
selbach, Y. Liu, and D. Mailly, Phys. Rev. Lett. 95, 
097004 (2005). 

V. O. Dolocan, P. Lejay, D. Mailly, and K. Hasselbach, 
Phys. Rev. B 74, 144505 (2006). 

N. A. Logoboy and E. B. Sonin, Phys. Rev. B 79, 094511 

(2009) . 

A. Bouhon and M. Sigrist, N. Jo. Phys. 12, 043031 

(2010) . 

F. Kidwingira, J. D. Strand, D. J. V. Harlingen, and 
Y. Maeno, Science 314, 1267 (2006). 



